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Abstract 

The hypothesis that the energy-momentum tensor of ordinary matter is 
not conserved separately, leads to a non-adiabatic expansion and, in many 
cases, to an Universe older than usual. This may provide a solution for 
the entropy and age problems of the Standard Cosmological Model. We 
consider two different theories of this type, and we perform a perturbative 
analysis, leading to analytical expressions for the evolution of gravitational 
waves, rotational modes and density perturbations. One of these theories 
exhibits satisfactory properties at this level, while the other one should be 
discarded. 



1 Introduction 

The Cosmological Standard Model is based on the General Relativity theory, cou- 
pled to a perfect fluid energy-momentum tensor, with geometry determined by 
the Robertson- Walker metric. It predicts an initial hot phase, when the baryoge- 
nesis and primordial nucleosynthesis take place, followed by the matter dominated 
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phase, when galactical structures are formed. This model can be considered as 
successful from many points of view, but it requires very stringent initial condi- 
tions. For example, the entropy observed today is of the order of 10 87 [|l|, which 
can be also expressed by saying that the actual density is very close to the critical 
density necessary for the existence of a spatially flat Universe. 

The thermodynamic equilibrium must also be imposed as an initial condition, 
at least while the nucleosynthesis takes place, due to the existence of the horizon 
distance smaller than the Hubble radius at that period. Moreover, the forma- 
tion of structures are well explained if we admit the mean density fluctuations 
of the order A oc 10~ 5 when the radiation decouples from matter; but there is 
no clear mechanism explaining the amplification up to this value of small initial 
perturbations, originating in the very beginning of the Universe. 

The inflationary model, based on the existence of a de Sitter phase prior to the 
nucleosynthesis, copes with the entropy and horizon problem. It also proposes a 
mechanism of generation and evolution of inhomogeneties. In what concerns the 
explanation of structure formation, inflation is in competition with other mecha- 
nisms, like those proposed by topological defect models. In contrast to the latter, 
it proposes gaussian perturbations of quantum mechanical origin. Until very re- 
cently, the observational results of the COBE satellite concerning the anisotropies 
of the Cosmic Microwave Background, seemed to favorize strongly the inflation- 
ary mechanism, but the gaussian nature of the perturbations has been questioned 
short time ago @. If this result is confirmed, it will become a serious drawback 
for the inflationary models. 

Various other cosmological models try to propose alternative solutions to the 
problems that remain open in the Standard Cosmological Model, without using 
the inflationary paradigm. Many of them introduce "exotic" types of matter 
sources, or modify the dynamical equations||. The models introducing a "cosmic 
string fluid" 0] can lead to good results concerning the horizon problem, but ex- 
hibit instabilities at perturbative level ||. Another interesting possibility consists 
in considering a "non conservative" theory of gravity, meaning a gravity theory 
in which the energy-momentum tensor does not conserve in the usual way. The 
propositions of this type go back to the work of Rastall 0, and have again at- 
tracted some attention recently ]7| ||. Such a scheme leads to an older Universe 
in comparison with the standard results, creating at the same time more entropy 
during its evolution. 

On a purely formal level, these "non-conservative" theories are equivalent to 
the introduction of exotic matter, whose main feature is negative pressure. More 
fundamentally, they lead to a radical modification of the dynamics of the Universe. 
We find frequently in the litterature statements explaining that such theories are 
incompatible with the equivalence principle. But the equivalence principle means 
that particles follow geodesic curves satisfying the equation u^u v .^ = 0, and this 
requirement is totally independent of T^ v '. ifM = 0, and there is no incompatibility 
between the assumption T^ u -^ ^ and the equivalence principle ||. 

Here we propose to study these theories from the perturbative point of view, 
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paying special attention to the evolution of gravitational waves, rotational modes 
and density perturbations. The perturbative analysis intends to give more infor- 
mation concerning the structure formation in the realm of these theories, checking 
their stability at the same time. 

In order to do this, we shall consider two different models. The first one, 
referred to as "Theory I", is based on the Brans-Dicke theory, rewritten in the 
Einstein frame. In this frame, matter is necessarily coupled to the scalar field. 
Certain considerations on the physical aspects of the two different frames have 
been presented in reference ||. In spite of serious objections against the theories 
of this kind, we will consider it here mainly as a reference model. The second 
model, referred to as "Theory II", is based on the original proposition due to 
Rastall, and in our opinion has better physical justification. The results obtained 
in the framework of this model indicate that the improvements with respect to the 
standard scenario are not compromised, at least at the perturbative level. This 
theory must be also implemented with a convenient mechanism for generation 
and proper evolution of inhomogeneties in an expanding globally homogenous 
and isotropic Universe. 

In the next section, we present briefly the two theories and discuss their cos- 
mological solution and some of its implications. In section 3, we perform a per- 
turbative study, considering perturbations of spin 2, 1 and 0. The comments and 
conclusions are contained in section 4. 

2 Generalized Einstein's theories 

The main point to be assured in order to solve the entropy problem without 
inflation, is to introduce a non adiabatic evolution of the Universe. This can 
be obtained by giving up the condition that the energy-momentum tensor has 
a vanishing covariant divergence. We can achieve this goal essentially in two 
different ways. First, we can start with a non- minimal coupling between the 
gravity and the scalar field on the one hand, and the matter on the other hand 
(like e.g. in the Brans-Dicke theory); then, through a conformal transformation 
we can recover the minimal coupling between the gravity and the scalar field, 
but with a more complicate coupling between matter, the metric field and the 
scalar field. The second possible choice is to modify Einstein's equations in such a 
way that the Bianchi identities are no longer satisfied in their normal form, with a 
conserved energy-momentum tensor on the right hand side, and in order to rectify 
the situation, we must modify the conservation laws for matter, which amounts 
to include new relationship between the gravitational and other forms of energy. 

We will call the first choice Theory I and the second one Theory II. We shall 
consider now separately each of these possibilities, analyzing their implications on 
the corresponding cosmological scenarios. 
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2.1 Theory I 

We start with the Brans-Dicke theory coupled to the matter in the ordinary way 

L=^/=gUR-u^-)+L m {^ . (1) 



Now, we redefine the field and perform a conformal transformation as in JlOj 



R-a, p aA+L m {e 2 ^g pp ) , (5) 



= t £ ~ 2 " a ' (2) 
9^ = e 20 °g^ , (3) 

= ■ (4) 

Combining the equations (0, ^ |]), we obtain a new lagrangian which can be 
written as : 

K 

where we have suppressed the tildes in order to simplify the notation. The varia- 
tional principle leads then to the following equations : 

1 1 

Rfiv - -j^g^R = kT im/ + a.^a-y - -g pu a- p a' p , (6) 

Ua = -KpT , (7) 
1^% = p<r"T , (8) 

where T = g^ u T^ v is the trace of the energy-momentum tensor. Imposing the 
Robertson- Walker metric 

ds 2 = dt 2 - ^-^-(dx 2 + dy 2 + dz 2 ) , (9) 
1 + f r 2 

and with the energy-momentum tensor of a perfect fluid, the above system reduces 
to the following three equations : 

3^ + 3^ = K P + \{a) 2 , (10) 
cr or 2 

a + 3-a = -K0(p-3p) , (11) 
a 

p + 3-(p + p) = P&(p-3p) . (12) 
a 

The solutions to these equations for the case k = are easily obtained : 
• a = —1: 

aocf , r=— , a = -—\nt ; (13) 
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a = 0: 



a = \: 



2 A3 

a(xf ' r= sTW ' a = -JTW lnt ; (14) 



Note that when — § < uj < oo, oo > /3 2 > 0. 

These solutions display the following characteristic features : 

1 . The solution for the radiative Universe is the same as in General Relativity, 
since the trace of the energy-momentum tensor is zero; 

2. The solution for a = — 1 exhibits an inflationary behaviour when uj > ~. 
The same situation occurs in the original Brans-Dicke theory [ |TT|j ; 

3. The solution for the material phase does not exhibit an inflationary behavior 
and the scale factor evolves always slower than in the standard model (0 < 
r < |). This implies that the age of the Universe in this case is less than in 
the standard case since t = rH , where H is Hubble's constant measured 
today 



Concerning the entropy production, using the equation (0), and the relation 
S = ^ p+ y " |J , where S and T are the total entropy and temperature of the 
Universe, we obtain the following solution for S, 



(3 2 

S = S exp da = S exp(- - - ^ In t) . (16) 

Hence the entropy always decreases with time during the material phase, while it 
remains constant during the radiative phase. 

2.2 Theory II 

The modified field equations are written with a a new parameter A as 

Rnv — -^9fj,vR = . (17) 



The divergence of (|T7|) leads to 

T^.,, = -^R U . (18) 

Obviously these equations coincide with the standard ones when A = 1. We can 
write this system of field equations in an equivalent form as 

1 / 7 — 1 \ 
Rtxv - ^g^R = k [T^ —g^T) , (19) 

T^ ]IX = l^-V" , (20) 
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where A = vp^-, and 7 = . 

One checks easily that the covariant divergence of the energy-momentum ten- 
sor vanishes and the equations become standard when 7=1 (i.e. with A = 1). 
This form of the system suggests that the modification amounts in fact to a re- 
defining the energy-momentum tensor by inclusion of part of gravitational effects 
directly, and does not affect the geometric framework of the theory. 

This type of modification is less radical than the departure from the Einstein- 
Hilbert lagrangian consisting in replacement of the scalar curvature R by a func- 



tion f(R), as in [12 1; the common feature remaining in both cases is that the extra 
geometrical terms may be interpreted as a modification of the energy-momentum 
tensor. 

It is generally argued that theories of this type can not be derived from a 
variational principle. However, we can define a matter lagrangian leading to the 
equations ([19], |20[). For example the matter lagrangian of the form : 

L m = \f II g(x(p+p)u ,J 'u u g fll/ + yp+ zp^j (21) 

leads to Einstein's equations when x — 1, y — — 1 and z = —3, and to the 
generalized Einstein's equations when we choose x — 1, y — —7 and z = —37. 

We look for solutions with spatially flat Robert son- Walker metric (H), where 
a(t) is the scale factor describing the evolution of the Universe. We consider the 
energy-momentum tensor of a perfect fluid, 

= (p + p)u»u v - pg» v , (22) 

where p, p and denote respectively, the energy density, the pressure and the 
four-velocity of the fluid. We assume a barotropic equation of state p = ap, 
where a shall take on the values —1, | and 0, corresponding to the vacuum 
energy density, radiation and matter dominated equation of state respectively, 
representing the most important cases for the primordial cosmological models. 

For these three equations of state we can easily determine simple solutions, 
namely, 

• a = —1: 



a 



1. 

3" 



a(t) oc expffi , H = \j^{3 - 2 7 p) , (23) 
p{t) oc constant ; (24) 



a(t) oc ti , p(t) oc r 2 ; (25) 



a = 0: 

a(t) oc t 1 "? , p{t) oc r 2 ; (26) 
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Among the above solutions only the case a = presents an improvement with 
respect to the standard model. The age of the Universe for this case can be 
evaluated as being equal to 

t = (1 - ^Ho- 1 , (27) 
leading to a Universe older than in standard mocdel for 7 < 1. On the other 



hand, using the equation (f20|), and the same relations between p, p, S and T as in 
the previous section, we can deduce the entropy growth in the matter dominated 
Universe according to the power law 

S = Sot 1 ^ . (28) 

Hence, the condition 7 < 1 implies at once an older Universe and an increase of 
the total entropy during the matter dominated phase. Moreover, the Universe 
exhibits an inflationary regime in the matter dominated era if 7 < 0. We can 
also note that due to the particular coupling between matter and gravity, it is 
possible to get negative values of the critical mass density parameter Q and at 
same time a closed or spatially flat Universe [0]. This feature is quite similar to 
the one displayed by a model which couples ordinary matter with a string induced 
fluid, when closed Universe can evolve dynamically as an open one (|J. However, 
it turns out that such models are hydrodynamically unstable ||. 



3 Perturbative analysis 

The Theory I considered above does not lead, in principle, to satisfactory results 
relative to the age of the Universe and the entropy generation. The Theory II, 
on the other hand, gives acceptable results if the parameter 7 is restricted to the 
values < 7 < 1. A perturbative analysis of these theories may give new insights 
into the physical consequences they lead to. In spite of previously shown negative 
results from the Theory I, we analyze it also from the perturbative point of view, 
for the sake of completeness. Our aim in both cases is to check if these theories 
can describe adequately the primordial gravitational waves, vectorial modes and 
density perturbations. In other words, we would like to see if they can lead to 
growing vectorial modes, associated with rotations, which is another drawback of 
the standard model. Also, if they can provide a mechanism for the amplification of 
density perturbations compatible with observations, since it is generally assumed 
that they are very small in the beginning of the Universe, and must be amplified 
during all its history in order to explain the observed large scale structures, as 
well as the observed distortions of the Cosmic Background Radiation. 

Technically, we use the standard procedure of Lifschitz-Khalatnikov [HJ. We 



isert into the field equations the perturbed quantities = g^ v + h^, p = p + 5p 
and p = p+8p, where g^, p and p are the background solutions around which small 
inhomogeneties are created. Using the invariance of field equations with respect 
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to the infinitesimal coordinate transformations, we may fix coordinate conditions 
eliminating unphysical degrees of freedom. Here we shall use the synchronous 
coordinate condition = 0. It is well known [R], that after that a residual 
coordinate freedom remains, leading to a non-physical degree of freedom in the 
density perturbations. In our final results, we will consider the physical modes 
only. The derivation of the perturbed field equations is standard. In addition to 
the considerations exposed above, we note that the perturbed quantities can 
be split into a sum of tensorial, vectorial and scalar modes as follows: 



h 



i.i 



h TT ij + h T m + h s m + Wfiij 



(29) 



where the first term is the transverse tracelless component (spin 2), the second 
one represents the vectorial mode (spin 1) and the last two ones are identified 
as scalar modes (spin 0). The derivatives are defined in the three-dimensional 
spatial section and 7^- is the induced metric of this spatial section. Moreover, it 
can be shown that the spatial components of the perturbed functions obey the 
eigenvalue equation V 2 Q — —k 2 Q, where k 2 = k 2 , and k is the wave vector of the 
perturbation which is supposed to propagate as a plane wave. 

We shall consider now the perturbed equations and their solutions in each of 
the aformentioned theories. In all the following formulae, k 2 = k 2 . 

3.1 Theory I 
• Spin 2: 



(30) 



Spin 1: 



Spin 0: 



'2 - 3a) - + /?<t (1 - 3a) 
a 



5u l = 



(31) 



h + 2-h-4fix-[H^) 



X 



a 2 )(l 



i + 3-rV^-/+3h)--^(l-3«)A 



.a 2 



3a)A 



A + (l + a)(6--/i)-/?(l-3a)x 



1 + a)Q + (1 + a) ( (2 - 3a)- + (3{l - 3a)&) 

+^ (aA +(3(1- 3a) x 









(32) 
(33) 
(34) 

(35) 
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In the expressions above, A = O = 8u\, h = ^ and \ = 5a. 

The equation describing the tensorial mode admits a simple solution when 
a = in terms of the conformal time defined as dt = adrj: 



oc Tj 



l + r 

2(l-r) 



c±J ±u (krj) 



(36) 



where v 



3-2f3 2 



>{l+2 ^ji J±v are Bessel functions with two possible values of index, 

±u, and c± are constants of integration. The results are essentially similar to 

i 



those obtained in the General Relativity framework for a — | [15], and in the 
original Brans-Dicke theory jTjJ for a = — 1. The vectorial mode equation also 
has a simple solution with 

(37) 



8u l oc t 3 + 2 <3 2 

Growing rotational modes can appear if j3 2 



— § < co> < — |, implying \ < r < 



> 1. This leads to the restriction 
~ - ,. ., ^ |. s It should be noted that such values of r 

correspond to a very young Universe. 

For the density perturbations (spin 0), the solution is exactly the same as 
for the radiative phase in General Relativity. As it is well known, there is no 
significant amplification of the perturbations during this phase fll5| . Nevertheless, 
unlike in General Relativity, we can have non-zero solutions for the perturbations 
during the de Sitter phase. Choosing a = — 1 in the perturbed equations, we find 
first that A = 4f3x- Using this relation, we can determine a third order diferential 
equation satisfied by A: 



+ 



,2- 



A+|52-?']A + 

a o / 



'-) + 3( - + (K) ~ -~)-2a 2 + 16/? 2 (3(^-) - -V ) A 



,d 2 



1 



a a 

-(-) 2 + 16/3 2 3-(2- - --) - -& 2 

a a \ a a a a a 



2 aa 



A = 



(38) 



As usual, A oc ]■ is a solution that is related to the residual coordinate freedom. 
Inserting A = l£ and using the conformal time such that dt = adr], then writing 
C = V p £) where p = |, we find the following equation for £: 



V 



1 - 



1,1 + 



4 V 1 - 8/3 2 ' r] 2 . 

This is a Bessel equation. Turning back to A, we have the solution: 



1 f 3 
A = J r 1 2 c±J±v(kr])dri , 



(39) 



(40) 



with v = ~ 1+ _ 8 ^ 2 2 and s 



8/3 2 -l • 
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This solution is not valid for n — 1 when this is the case, we find an Euler 
equation,: 

A+6y + (A; 2 + 6)^ + A; 2 A = . (41) 
Discarding the solution connected with the residual coordinate freedom we have : 

A = r i±>/i=fc* . (42) 

The main feature of these solution is that asymptotically, for small k (large wave- 
length of the perturbations), it does not contain growing modes. 

For a = 0, we can find solutions in the long wavelength limit, k — > 0. They 

Aoct " • p ^(- i+w!±(5 - 2fl ) • <43) 

The maximum value of p appears when to — > oo, leading to p — |. In this limit 
the theory coincides with General Relativity. Hence, during the material phase, 
the perturbations grow less rapidly than in the standard model. 

3.2 Theory II 

For this theory, the perturbed equations read 
• Spin 2: 



• Spin 1: 



Spin 0: 



h -%+((-) -2^)h l3 = ; (44) 

(X \ (Jj Ub / 



(1 + p)5u i + (1 + (3){2 - 3/3) -Su* = ; (45) 

CL 



h + 2-h = 6(^)(1 + /3)A , (46) 

A + (l + P)(Q-h) = , (47) 

(1 + ^)6 + (1 + ^(2-3^-0 = /3(^)A , (48) 

(x a 



where here (3 = ^E^z^z^ ■ Hence, the equations (so, the solutions) are similar 
to those of the standard model, but with a "modified equation of state" in which 
(3 replaces a. We can verify that, for a = — 1, (3 = —1, a — |, (3 — |, and 
a — 0, f3 — |+^. The results are the same as in the standard model for the 
false vacuum and radiative phases. During these periods, the rotational modes 
decrease (or simply does not exist, as in for the false vacuum case), and the density 
perturbations exist only in the radiative case, but are amplified very moderately. 
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Let us analyze now the case a = 0. For < 7 < 3 (so that the weak energy- 
condition is satisfied), — | < (3 < 00. For < 7 < 1 (growing entropy and older 
Universe), we have — | < (3 < (the effective pressure is negative, but the strong 
energy condition is not violated). We consider now the solutions for the three 

different modes separately. The scale factor behaves as a oc or, in terms of 

3-7 

the conformal time, a oc 77 7 , 

• Spin 2: The solutions can be written as, 

1 6 — 1 6 — 37 

oc 772 7 J u (krj) , v = — . (49) 

27 

For large wavelengths, it behaves asymptotically as 

6-27 

hij oc tj 7 . (50) 

In the interval of interest, the gravitational waves are amplified, but not in 
a exponential way. 

• There are growing rotational perturbations during the material era provided 
that 

3/?-2>0^/?>^ 7 > j| • (51) 
However, this leads to an young Universe and a decreasing entropy(7 > 1); 

• For the density perturbations, we can note that: 

1. For 7 > 1, the solutions are equivalent to those of the standard model 
with an equation of state p = ap, a > 0, namely, 

Aoc — / v -2 C± J ±u {kr]) , v = . (52) 

777 J ^7 

2. For < 7 < 1 and 7 < 0, the solution are the same as in General 
Relativity, with — | < a < 0, which amounts to replacing the Bessel's 
functions by the modified Bessel's functions (see, e.g. JH]]) : 

A oc 4r / V 1 (cilu(k v ) + c 2 K u (k V )) . (53) 
fj 7 J \ / 

In both cases, the long wavelength limit yields the same asymptotic be- 
haviour 

A oc rf oc th . (54) 

We can see that the perturbations, in this limit, and with 7 > 0, grow faster 
than in the Standard Model when the Universe is younger, and slower when 
the Universe is older. In both cases, we can expect the final amplification of 
the perturbations to be of the same order. If 7 < 0, the density perturbations 
are decreasing. 



11 



4 Conclusions 



Models where the momentum-energy tensor is not conserved per se may give rea- 
sonable results while coping with some traditional problems posed by the Standard 
Cosmological Model. Here we have considered two kinds of such theories. The first 
one is based on a non-trivial interaction between the matter and the scalar field, 
which can be obtained from the Brans-Dicke theory through conformal transfor- 
mation. The second one is based on a direct modification of Einstein's theory, via 
introduction of the parameter 7 in Einstein's equations. 

The first theory is marred by problems from the Cosmological point of view 
even in the background solution, since it predicts an Universe younger than the 
one predicted by the standard model, and decreasing entropy during the matter 
dominated phase. The second theory, in contrast, gives the opposite result when 
< 7 < 1. For 7 < it predicts an inflationary behaviour during the matter- 
dominated phase. 

Turning to the structure formation problem, we can state the following results. 

1. In the realm of the first theory, the gravitational mode has qualitatively the 
same behaviour than in the standard model. The rotational mode exhibit 
a growing solution but for a very young Universe. But, contrary to the 
standard model, there exist density perturbations in the false vacuum case, 
which, however, have a decreasing behaviour in the long wavelength limit. In 
the radiative phase, the behaviour is the same than in the standard model, 
but during the matter dominated phase, the density perturbations grow less 
rapidly, which is even more serious drawback if we remember that the age 
of the Universe is also lower than in the standard case; 

2. In the second theory, we have the same behaviour that within the standard 
model if < 7 < 1 for gravitational and rotational modes. Rotational modes 
may exhibit a growing tendency for 7 > | but this leads to a very young 
Universe and the entropy decreasing during the matter dominated phase. 
The behaviour of the density perturbations remains unchanged with respect 
to the standard results during the false vacuum and radiative phases. But 
during the matter dominated phase there is a fast growing mode (in the long- 
wavelength limit) if 7 > 1, and a slower one if 7 < 1. The amplification is 
essentially the same in this phase due to the modifications in the evaluation 
of the age of the Universe. 

The first theory does not lead to acceptable results neither for the background 
solution, nor at the perturbative level, whereas the second one gives interesting 
background solutions and, at least, is not disqualified by the perturbation analy- 
sis. This can be judged rather encouraging, since as in the standard model, the 
second theory needs a mechanism for generation and amplification of perturba- 
tions in phases prior to the matter dominated era; but, on the other hand, we 
remember that certain models which give interesting background solutions are 
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entirely compromised at the perturbative level, as is the case of the models based 
on string like fluid ||, and which is not the case here. 

For 7 < 1, this theory leads to a reasonable entropy production, to an older 
Universe and to growing modes for density perturbations in the material phase. 

Is it possible to extract some observational limit for the parameter 7 from the 
above perturbative analysis? In principle, yes. We could employ for example the 
observations on the anisotropies of the cosmic microwave background, measured 
by the COBE satellite. However, in order to do that, we must have a complete 
scenario at our disposal, including the above mentioned mechanism for the origin 
and amplification of the perturbations. In our case, this should be done without 
the help of inflation in order to be coherent with the proposals of the model. 

Even without that, we can extract some observational informations. One im- 
portant observational parameter is the spectral index n, characterizing the spec- 
trum of the perturbations. If n = 1 we have the so called Harrison-Zeldovich spec- 
trum: all perturbations, disregarding their wavelength, have the same amplitude 
when they enter the horizon. The COBE's results indicate a nearly Harrison- 
Zeldovich spectrum, with n = 1 representing possible value. 

The spectral index can be evaluated from the two point correlation function. 



It gives the following result [19, 18 



dlnS 2 

" = 1 + ^ • (55) 

where 5 2 is proportional to the two point correlation function, 

\A(k,t)\^k s J-^^(A(x,t)A(0,t)) . (56) 



Using the solutions fl52|), considering the long wavelength limit with the solution 
that leads to the Bunch-Davies vacuum (for details, see |l9j), we have 

5 oc k 3 ~ 2v , (57) 

leading to 

n = l + 6 ilzll . (58) 

7 

The observational limits of n seem to indicate a value very close to n = 1. Taking 
for simplicity (since we are only making a crude estimation) 0.8 < n < 1.2, we 
obtain 0.96 < 7 < 1.03. Following this simplified analysis, we find that the 
Universe today must be, from the Theory II, a little older than the one predicted 
by the standard model, with a very low production of entropy during the matter 
dominated phase. 
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